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Abstract 



The vacuum oscillation solution to the solar neutrino problem predicts char- 
acteristics variations of the observable neutrinos rates, as a result of the L/Ey 
dependence of the Vg. survival probability (L and E^, being the neutrino path- 
length and energy, respectively). The -Ej^-dependence can be studied through 
distortions of the recoil electron spectrum in the SuperKamiokande experi- 
ment. The -L-dependence can be investigated through a Fourier analysis of 
the signal in the SuperKamiokande and Borexino experiments. We discuss in 
detail the interplay among such observable variations of the signal, and show 
how they can help to test and constrain the vacuum oscillation solution(s). 
The analysis includes the 374-day SuperKamiokande data. 



PACS number(s): 26.65.+t, 13.15.-Fg, 14.60.Pq 
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I. INTRODUCTION 



Neutrino flavor oscillations [|I| with wavelength comparable to the Earth-Sun distance 
represent a solution to the deficit of solar z/'s 0] observed in the four pioneering 
underground experiments Homestake 0, Kamiokande 0, SAGE 0, and GALLEX p|, as 
compared to the standard solar model predictions [Q. The recent SuperKamiokande data 
p!0| confirm the deficit, and can be interpreted within the vacuum oscillation hypothesis 



as well |Tl|JT^ . The planned Borexino solar neutrino experiment |I3|JT^ (in construction), 
designed to detect monochromatic ''Be neutrinos {E^, = 0.86 MeV), is expected to test this 
hypothesis with unprecedented sensitivity ||T^|-|T^. 

An update of the current neutrino fiux measurements P, |T0| , p!8| -^ is given in Table I. 
Figure 1 shows our vacuum oscillation fit to the data of Table I, as obtained from a 
analysis (including solar model uncertainties as in We have assumed, for simplicity, 

two neutrino families. It can be seen that four regions (A, B, C, and D) are allowed at 
95% C.L., the absolute minimum being located within the solution B (Xmin — 3-4 and 
A^DF = 3 = 5-2). 

Various tests can be envisaged to discriminate among the four solutions in Fig. 1. In 
particular, since the neutrino oscillation length is proportional to the pathlength-to-energy 
ratio L/Ey, deviations of event distributions from the expected shape in either L or Ej^ (or 



related parameters) represent direct tests of neutrino vacuum oscillations (see, e.g., p2|-p^). 
It turns out that, in general, SuperKamiokande is more (less) sensitive than Borexino to 
£^j,-related (L-related) spectral shape deviations; therefore, the two experiments provide 
complementary tools to study the vacuum oscillation hypothesis. The purpose of this work 
is to investigate in detail the tests of the vacuum oscillation hypothesis that can be performed 
at SuperKamiokande and Borexino, and their interplay. The results will be shown in a form 
that makes easy to derive the experimental accuracy needed to perform a specific test. 

The plan of this paper is as follows. In Sec. II we discuss the tests of energy spectra 
deviations, in the light of the recent SuperKamiokande results. In Sec. Ill we discuss in 



detail the Fourier analysis of the signal, building upon our previous work [25|. In Sec. Ill 



we apply these tests to SuperKamiokande and Borexino, both separately and jointly. In 
Sec. IV we draw our conclusions. Some technical aspects of our analysis are described in 
the Appendix. 



II. £;^-RELATED TESTS: AVERAGE ELECTRON KINETIC ENERGY 

Both SuperKamiokande and Borexino can measure the energy spectrum of recoil elec- 
trons from neutrino scattering. The standard (i.e., no oscillation) SuperKamiokande spec- 
trum can be found in Fig. 4 of Ref. [|ll|]. For completeness, we show in Fig. 2 the standard 
electron energy spectrum in Borexino (details about the inclusion of energy threshold and 
resolution effects are given in the Appendix). The main contribution to the spectrum in 
Fig. 2 is given by the 0.86 MeV ^Be line, which is responsible for the Compton edge at 
~ 0.66 MeV. The edge is smeared by the finite energy resolution. The rise at low energies is 
due to pp neutrinos, while the tail at high energies is basically due to CNO neutrinos. The 
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prospective analysis window is also indicated in Fig. 2. The standard neutrino fluxes have 
been taken from 

Distortions of the neutrino energy spectrum due to oscillations are generally reflected (al- 
though somewhat degraded) in the electron energy spectrum. An effective parametrization 
of such distortions is given by the fractional variation of the average kinetic energy (T) of 



the electron, an approach extensively developed in p2[ and applied to the SuperKamiokande 
In particular, the analysis of the most recent (374 day) measurements 



data in ITT 



of the electron spectrum in SuperKamiokande [p!0| ,p^^-|3T| gives a fractional deviation 



A(T) 
(T) 



X 100 = 0.95 ±0.73 



(1) 



Such deviation is consistent with zero (no oscillation) at the level of 1.2a, and disfavors 
scenarios predicting negative values for A(T)/ (T). 

Concerning Borexino, the expected shape variations of the energy distribution are very 
small, since the main contribution to the electron spectrum comes from a monoenergetic 
source of neutrinos (^Be), rather than from a continuous source as in SuperKamiokande 
(^B). Neutrino oscillations are expected to change significantly the global Borexino rate but 
not its energy spectrum. 



III. L-RELATED TESTS: FOURIER EXPANSION OF THE SIGNAL 

The variations of the solar neutrino pathlength L due to the eccentricity of the Earth 
orbit produce a geometrical {l/L"^) modulation of the neutrino flux. Additional semiannual 



modulations are expected in the presence of neutrino oscillations |3^. The Fourier analysis 



of the measured flux represents an effective tool to study both kinds of L-related modulations 
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In this Section we outline the Fourier analysis of the signal observable in solar neutrino 
experiments. In the first three subsections, we describe the general formalism and the results 
for the "no oscillation" and "2z/ oscillation" cases (see also Ref. [^] for further details). In 



the last two subsections we generalize the analysis to Su oscillations and then discuss a useful 
check of both the symmetry properties and the estimated uncertainties of the signal. 



A. General formalism 

The Earth orbit radius, L, varies periodically in time (t) around its average value, Lq = 
1.496 X 10® km, according to 

L(t) = Lo(l-ecos^)+0(e2), (2) 

where e = 0.0167 is the orbit eccentricity and T = 1 yr (t = at the perihelion). Terms of 
0(£:^) or higher are negligible for our purposes. 

The neutrino signal S is also, in general, a periodic function. Assuming a constant 
background B, the total observed neutrino rate R is: 
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R{t) = B + S{t) 



(3) 



For symmetry reasons, the analysis can be restricted to the time interval [0,T/2].Q It is 
understood that events collected in subsequent half-years must be symmetrically folded 
in this interval. The data sample consists then of events collected at different times 
{ti}i<i<N, with ti G [0, T/2] and N equal to the total sum of background and signal events, 
= Ns + Nb- Notice that, in general, one can determine Nb and Ns but cannot distinguish 
background and signal on an event-by-event basis. 

The expansion of the signal in terms of Fourier components reads 

/ °° 27Tr}f\ 

S{t) = s[l + 2^Ucos'-^j , (4) 

where S is the time-averaged signal 

2 /-^/s 

S = -J^ dtS{t). (5) 

The n-th harmonic corresponds to a period of l/n yr. The explicit form of /„ reads 

2 2iTnt 

fn=gfj^ dtR{t) COS ^ (6) 

= jr cos— ^ (0<t.<T/2), (7) 

where Eqs. (|]) and (0) represent the theoretical definition and the experimental determina- 
tion of the /n's, respectively p5 |. 

Assuming purely statistical fluctuations of the signal and of the background, the variance 
of the /„'s reads: 

... 1 + f2n + Nb/Ns 

var(/„) = . (8) 

It turns out that the values of the /2n's are <C 1 in all cases of practical interest. Therefore, 



to a good approximation, the one-sigma statistical error o"/ = y var(/„) affecting /„ is given 
by 



for any n. The correlations between the statistical errors of different harmonics are also 
negligible ||2^ . 

Finally, the general expression of the signal S expected in the presence of oscillations is 



Experimental tests of the symmetry properties of the signal are discussed at the end of this 
section. 
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^^JdE XiE) [a,{E)P{E, t) + a^{E){l - P{E, t))] , (10) 

where E is the neutrino energy, A is the neutrino energy spectrum, cTg (a^) is the z/g (z/j., 
X = n, t) interaction cross section, and P is the z/g survival probabihty, which varies in time 
through L(t) ||32|. It is understood that the cross sections dg^x must be corrected for energy 
threshold and resolution effects, as described in the Appendix. 

B. Standard (no oscillation) case 

In the standard (no oscillation) case, characterized by P{E,t) = 1, the signal S varies 
as S'^^'^{t) oc L~'^{t). The standard Fourier components are simply given by 

= ^5nl , (11) 

i.e., only the first harmonic is nonzero and measures the Earth's orbit eccentricity. 

C. 2u oscillation case 

In this case, Uf, is a linear combination of two mass eigenstates (z/i, z/2) characterized by 
a mixing angle 6 

Vf. = cos 61/1 + sin 6 z/2 (12) 
and by a squared mass difference (5m^, 

= \ml — mil . (13) 
The corresponding z/g survival probability is given by 

P-(E, t) = l-l sin^ 26 (1 - cos ^!^) . (14) 



The Fourier components can be cast in the following, compact form pS 



1 - sin^ 2eDo{6m'^j 
where the (detector-dependent) functions Dn are given by 



r2u _ g^^ni - sin^ 26 Dn{6m^) 



/ dE A (cTg — Gx) Un 

Dn{6m^) = ^ J {n > 0) (16) 

2 dEXa, 



and the universal (i.e., detector-independent) functions [/„ are given by 



Un{z) = COS (^Z 



JJez) 



T 

where z = 5m?LQ/2E and J„ is the Bessel function of order n. Notice that, although our 
calculations are of 0(e), all orders in ez are kept, since z may be large. 
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D. 3z/ oscillation case 



In this case, Ve is a linear combination of three mass eigenstates (i^i, 1^2, t's), usually 
parametrized in one of the following two forms: 



(17) 
(18) 



where s = sin, c = cos, and U^^ + f^e2 + ^( 
hierarchy of mass differences: 



2 

e3 



1 (see, e.g., |2^). In addition, we assume a 



brv? <^ — 2I 



(19) 



The hypothesis ([19|) covers most of the situations of phenomenological interest . In this 
case, the 3z/ oscillation probability is given by 



P^^E, t) = 1 - 2 Ull\ - Ul^) - 2 Ul^Ul^ (1 - cos 



6m^L{t)' 



2E 



We have worked out the corresponding Fourier components, which read: 



J n 



1 - 4 [7^3(1 - f/23) r - 4 Ul,U% Do{6m^ 



1 - 44>r - sLclDo{6m^) 
where the functions D„ are defined in Eq. (|16D and 

J dE X{(7e — (Jx 



2 / dE \a. 



(20) 
(21) 

(22) 
(23) 



(24) 



Notice that f^" f^" for U, 



0, as it should. 



E. A useful consistency check 

A priori, the signal S{t) must obey the symmetry S{t) = S(T — t) in the one-year 
interval [0, T] (either with or without oscillations). We have made use of this property in 
Eqs. (1^-0). The experimental test of such symmetry property is not without merit, since 
its failure might signal systematic, time-dependent effects, such as unexpected variations of 
the detection efficiency or of the background level. Within our approach, this is equivalent 
to check that the "sine" Fourier components Qn, defined as 

9n = j;r sin— ^ (0<t.<T), (25) 
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are identically zero, as they should (both in the standard case and in the presence of possible 
oscillations): 



qT = 9T = . (26) 

Notice that, in Eq. (pSj) , the event "arrival times" are folded in the interval [0, T] and not 
in [0, r/2]. 

The experimentally inferred will be distributed around zero with a variance var(5f„). 
As far as statistical fluctuations are concerned, the calculation of var(5f„) is analogous to 



var(/„,) and gives: 



var(^„) = . (27) 

However, as already noticed, the /2n's are generally ^ 1 and thus can be neglected in the 



above equation. Therefore, the statistical uncertainty (jg = y var((7„) is approximately equal 
for all the gnS and has the same expression as o"/ [Eq. (|^)]. 

In conclusion, it is useful to check that the semiannual variations of the solar neutrino 
signal are indeed symmetric in time. This implies that the "sine" Fourier components gn 
defined in Eq. (^) should form a distribution with mean value (gn) = and standard 

deviation ag = ^ (Ns + Nb) /2A^|. This test is universal, i.e. it is valid with or without 
oscillation effects. Any deviation of the mean and variance of the gn distribution from their 
standard values (0 and cr^, respectively) would indicate the presence (and the magnitude) of 
systematic effects beyond the purely statistical fluctuations of either the background or the 
signal. These effects (if any) should then be accounted for by the experimentalists, before 
performing an unbiased analysis of the time variations of the signal. 



IV. VACUUM OSCILLATION TESTS AT SUPERKAMIOKANDE AND 

BOREXINO 



In this section we analyze the tests of the vacuum oscillation hypothesis that can be 
performed at SuperKamiokande and Borexino, both separately and jointly. We consider 
two observables for each experiment. In particular, we analyze A(T)/(T) and /i — £ for 
SuperKamiokande (SK), and fi — e and /2 for Borexino (BX): 

Observables | ^^P^^^^ 7%'''^' " (SuperKamiokande) , 

I ff^ — e and (Borexmo) . ^ ^ 

All the above variables are zero in the standard (no oscillation) case. In the oscillation range 
of interest, Fourier components with n > 1 {n > 2) are not relevant for SuperKamiokande 
(Borexino) For simplicity, we will consider only 2z/ oscillations, and the corresponding 
preferred regions A, B, C, and D of Fig. 1. 
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A. Tests at SuperKamiokande 



Figure 3 shows the four solutions A, B, C, and D (gray regions) and the no oscillation 
point (star) in the plane charted by the parameters A(T)/(T) (the fractional deviation of 
the mean electron kinetic energy) and /f ^ — e (the deviation of the first Fourier component 
from its standard value). Also shown is the horizontal band allowed at ilo" by the A(T) / (T) 
datum of Eq. (|l]). 

Figure 3 evidences that solutions C and D (which predict large, negative values for 
A(T)/(T)) are highly disfavored by the SuperKamiokande measurement of Eq. ([l|). In 
particular, solution C is disfavored at > 4(T and solution D at > 6cr. On the other hand, the 
datum of Eq. (0) is unable to discriminate among the solutions A, B, and the no oscillation 
point at the 2a level, although solution A seems to be preferred. We can summarize these 
findings by saying that, under the hypothesis of 2v oscillations, the combined data of Table I 
and Eq. (|I|) select the solutions A and B in Fig. 1, corresponding to the following approximate 
ranges (at 95% CL.) for the neutrino mass-mixing parameters: 

0.59 < < 0.84 (xlO'io eV^) , (29) 
0.66 < sin^ 26 < I . (30) 

Notice that the spread of the above parameters is only about ±20%. Also notice that 2i/ 
maximal mixing (sin^ 26 = 1) is allowed only in a restricted range of (5m^ (~ 0.59-0.61 x 10~^° 
eV2). 

Concerning /i, we cannot infer its value from the limited data which are publicly available 
T0| , |27| -pl| . At any rate, the estimated uncertainty aj of /i [see Eq. (||)] for ~ 1 yr of data 



taking is about ±0.02 [|TT| — too large to discriminate any of the solutions A, B, C, and D. 
Significantly higher statistics are needed to reduce such error. 

An interesting feature of Fig. 3 is the tight correlation between the variables fi — e and 
A(T)/(T), which parametrize L-related and i?t,-related spectral distortions, respectively. 
Such correlation can be traced to the L/Ei, dependence of the oscillation probability, as 



emphasized in One can use such correlation to "predict" the value of /i — e in Su- 

perKamiokande for a given value of A(T)/(T) in Eq. ([^). More precisely, from Fig. 3 one 
derives that the values of fi — e compatible with both the horizontal band and the solution 
A should lie in the range ~ 0.005-0.01. In other words, one expects an enhancement of the 
semiannual modulations of the neutrino flux (relatively to the purely geometrical one) in the 
range [0.005/e, 0.01/e] = 30-60% at ~ la. This is a clear prediction that needs, however, 
several years of data taking at SuperKamiokande to be tested. 



B. Tests at Borexino 



Figure 4 shows the four solutions A, B, C, and D (gray regions) in the plane charted by 
the parameters /f ^ — e and (the deviation of the first two Fourier components from 
their standard values). The no oscillation point (star) corresponds to the origin. 

The amplitude of the second harmonic appears to be generally smaller than the first; 
nevertheless, both should be detectable in a sample of a few thousand events. E.g., for 
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Ns — Nb — 5000 events, the expected statistical error of /i and /2 is only about ±0.014 
[see Eq. (^], ensuring clear detection of semiannual modulations (provided that systematics 
do not dominate the error budget). However, the four solutions A, B, C, D are rather close 
to each other in the Fourier parameter space, and it might be difficult to distinguish among 
them using only these two variables. 

By comparing Figs. 3 and 4, it can be noticed that the solutions A and B (currently 
favored by SuperKamiokande data) predict rather different values for the first Fourier har- 
monic in SuperKamiokande and Borexino, as a result of the different energy ranges probed 
by these two experiments. In particular, within solution A it is always /f^ — e > 0, while 
ff^ — 6 can be either positive or negative. Therefore, at present the sign of the semiannual 
variations in Borexino is unpredictable. 

It is interesting to notice that, both in 2z/ and 3i> oscillations, the Fourier component 
ratio /2/(/i — s) depends only on Sm'^ and not on the mixing angle(s) [see Eqs. ([I5| ) and 
(Bl)]. Therefore, such ratio can constrain the value of in a model-independent way. 



C. Combination of SuperKamiokande and Borexino tests 

Figure 5 shows the four solutions A, B, C, and D (gray regions) in the four planes 
charted by the Borexino observables ff^ — e and f^^ (x-axes) and by the SuperKamiokande 
observables /f^ — e and A(T)/ (T) (y-axes). Also shown is the horizontal band allowed at 
±lcr by the A(T)/ (T) datum of Eq. (|l|). Within such band, the values of ff^ — e (as well 
as those of f2^) can be either positive or negative, as observed in the previous subsection. 
Figure 5 shows that the present indeterminacy in the sign of ff^ — e cannot be resolved by 
increasing the accuracy of the SuperKamiokande data. 

The similarity between the upper and lower panels in Fig. 5 is due to the tight correlation 
between the SuperKamiokande variables /f^ — e and A(T)/(T). This similarity should be 
reflected in the experimental data, if vacuum oscillations indeed occur. When experimental 
data will be available for all the four observables charting the panels in Fig. 5, one of the four 
solutions should be easily spotted in at least one of the four panels. However, in the unlucky 
situation of data points close to the no oscillation case, it might be difficult to distinguish 
between such case and solution B. The possibility to separate the no oscillation point from 
solution B would then depend decisively on the reduction of the A(T)/ (T) uncertainty (see 
upper panels of Fig. 5). In any case, the reader can judge the discriminating power of the 
two experiments by drawing prospective data points and error bars in each panel of Fig. 5. 



V. SUMMARY AND CONCLUSIONS 

We have studied tests of vacuum oscillations of solar neutrinos in the SuperKamiokande 
and Borexino experiments, both separately and jointly. The tests are sensitive to either 
energy or time variations of the neutrino flux. The interplay between such tests has been 
investigated, in the light of the most recent data (374 day) from the SuperKamiokande 
experiment. The results have been displayed in a graphical form (Figs. 3-5) that allows 
to determine easily the experimental accuracy needed to test the vacuum oscillation solu- 
tion(s). We have found that: (i) The total neutrino rates measured by solar v experiments 
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can be fitted in four distinct regions of the mass-mixing parameter space; (ii) Two of the 
four solutions are strongly disfavored by the SuperKamiokande energy spectrum data; [iii) 
The energy spectrum data do not discriminate significantly (at present) the remaining two 
solutions between them and from the no oscillation case; {iv) The amplitude of semian- 
nual variations of the solar u flux in SuperKamiokande is predicted to be about 30-60% 
in excess of the purely geometric one (at la); [v) The sign of semiannual variations (due 
to oscillations) in Borexino is not determined by present data; [vi) The joint information 
coming from the energy spectrum data (SuperKamiokande) and from the Fourier transform 
of the solar neutrino rates (SuperKamiokande, Borexino) can provide powerful tests of the 
vacuum oscillation solutions. Besides, we have generalized the Fourier transform formalism 
to three-flavor oscillations, and we have discussed a useful check of the time symmetry of 
the signal. 
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APPENDIX: ENERGY THRESHOLD AND RESOLUTION EFFECTS 



In the calculation of the expected signal [Eq. (p!0D ], it is understood that the Ua-e cross 
sections aa{E) (a = e, x) have to be properly corrected to take into account the detector 
energy resolution and the analysis window for each experiment. Here we give some details 
about such corrections. 

Both in Borexino and in SuperKamiokande, the finite energy resolution due to the photon 
statistics implies that the measured kinetic energy T of the scattered electron is distributed 
around the true kinetic energy T' according to a resolution function R{T, T') of the form 



R{T, T') 



1 



exp 



(T-T 
2^ 



r\2 



(Al) 



where 



s = soJT'/MeV 



(A2) 



and So = 57.7 KeV and 0.47 MeV for Borexino [|l^ and SuperKamiokande |T^,|TT]], respec- 
tively. On the other hand, the distribution of the true kinetic energy T' for an interacting 
neutrino of energy Ei, is dictated by the differential cross section daa^Ei,, T')/dT', that we 
take from 34 . The kinematic limits are: 



< T' < T'{E,) , T'{E, 



1 + me/2E, 



(A3) 
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Concerning the measured kinetic energy, the present analysis window for Su- 
per Kamiokande is [Tmin, ^max] = [6-5 McV — 20 McV — me]. The prospective analysis 
window for Borexino (as used in this paper) is [0.25, 0.8] MeV. 

For assigned values of Sq, T^in, and Tmax, the corrected cross section aa{Ey) is defined 

as: 



where the function W{T'), which embeds the detector specifications Sq, Tmin, and Tmax, is 
given by 




(A4) 



It is useful to reorder the integrands, obtaining 




(A5) 




(A6) 



with s as in Eq. 



and 




(A7) 
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TABLES 



TABLE L Neutrino event rates measured by solar neutrino experiments, and corresponding 
predictions from the standard solar model The quoted errors are at la. 



Experiment 



Ref. 



Data ib(stat.) ib(syst. 



Theory g 
9.3+^-^ 
6.62 
137 
137 
6.62 



Units 



Homestake 
Kamiokande 
SAGE 
GALLEX 
SuperKam. 



1 

i 

121 

m 



2.56 ±0.16 ±0.15 
2.80 ±0.19 ±0.33 
69.9l?:fti? 
76.4 ± 6.3l^-^ 

h0.06+0.d9 



2.37 



-0.05-0.06 



1.4 
-0.93 
-1.12 
+8 
7 
+8 
-7 
0.93 
-1.12 



SNU 
10^ cm-^g-i 

SNU 

SNU 
10^ cm-^s"^ 



12 



REFERENCES 



[1] B. Pontecorvo, Zh. Eksp. Teor. Fiz. 53, 1717 (1967) [Sov. Phys. JETP 26, 984 (1968)]; 

V. Gribov and B. Pontecorvo, Phys. Lett. 28B, 493 (1969). 
[2] S. L. Glashow and L. M. Krauss, Phys. Lett. B 190, 199 (1987). 
[3] P. L Krastev and S. T. Petcov, Phys. Rev. D 53, 1665 (1996). 

[4] J. N. Bahcall, Neutrino Astrophysics (Cambridge University Press, Cambridge, England, 
1989). 

[5] R. Davis (Homestake Collaboration), Prog. Part. Nucl. Phys. 32, 13 (1994). 

[6] Y. Fukuda et al. (Kamiokande Collaboration), Phys. Rev. Lett. 77, 1683 (1996). 

[7] J. N. Abdurashitov et al. (SAGE Collaboration), Phys. Rev. Lett. 77, 4708 (1996). 

[8] W. Hampel et al. (GALLEX Collaboration), 288, 384 (1996). 

[9] J. N. Bahcall and M. H. Pinsonneault, Rev. Mod. Phys. 67, 781 (1995). 
[10] Y. Totsuka (SuperKamiokande Collaboration), in the Proceedings of the 35th Interna- 
tional School of Subnuclear Physics (Erice, Italy, 1997), to appear. 
[11] G. L. Fogh, E. Lisi and D. Montanino, University of Bari Report No. BARI-TH/284-97, 

|hep-ph/970947^ , to appear in Astropart. Phys. 
[12] N. Hata and P. Langacker, Phys. Rev. D 56, 6107 (1997). 

[13] Borexino Collaboration, C. Arpesella et al.., "INFN Borexino proposal," Vols. I and II, 
edited by G. Bellini, R. Raghavan, et al. (University of Milan, 1992); J. Benziger, F. 
P. Calaprice et al., "Proposal for Participation in the Borexino Solar Neutrino Experi- 
ment," (Princeton University, 1996). 
[14] Updated information about the Borexino experiment can be found at the URL 
|http:/ / almime.mi.infn.it 



[15] Borexino Collaboration, G. Bellini et al. in TAUP '95, Proceedings of the 4th Interna- 
tional Workshop on Theoretical and Phenomenological Aspects of Underground Physics, 
Toledo, Spain, edited by A. Morales, J. Morales, and J. A. Villar [Nucl. Phys. (Proc. 
Suppl.) B48, 547 (1996)], p. 363; G. Bellini, in La Thmle '97, Proceedings of the 11th 
Rencontres de Physique de la Vallee d'Aoste, La Thuile, Italy, 1997, ed. by M. Greco 
(Frascati Physics Series, Vol. 9), p. 27. 

[16] R. S. Raghavan, Science 267, 45 (1995). 

[17] J. N. Bahcall and P. I. Krastev, Phys. Rev. C 55, 929 (1997). 

[18] K. Lande (Homestake Collaboration), in u^. Proceedings of the 4th International So- 
lar Neutrino Conference, Heidelberg, Germany, 1997, ed by W. Hampel (Max Planck 
Institute, Heidelberg, 1998), p. 85. 

[19] V. N. Gavrin (SAGE Collaboration), in News about SNU's, International Workshop on 
Solar Neutrinos (Institute of Theoretical Physics, Santa Barbara, CA, 1997). Proceed- 
ings available at the URL |http : / / www. itp . ucsb . edu / online / snu 



[20] W. Hampel (GALLEX Collaboration), in News about SNU's g 
[21] G. L. Fogh and E. Lisi, Astropart. Phys. 3, 185 (1995). 
[22] J. N. Bahcall, P. I. Krastev, and E. Lisi, Phys. Rev. C 55, 494 (1997). 
[23] B. Paid, G. L. Fogh, E. Lisi and D. Montanino, Phys. Rev. D 55, 1353 (1997). 
[24] S. P. Mikheyev and A. Yu. Smirnov, ICTP Report No. IC-97-113, |hep-ph/ 9708403 
[25] G. L. Fogh, E. Lisi and D. Montanino, Phys. Rev. D 56, 4374 (1997). 



13 



[26] G. L. Fogli, E. Lisi and D. Montanino, University of Bari Report No. BARI-TH/296-98, 
|hep-ph/9803309| , submitted to Phys. Lett. B. 

[27] M. Nakahata, in P4 '97, Proceedings of the APCPT Pacific Particle Physics Phe- 
nomenology Conference, (Seoul, Korea, 1997); also available at the URL [http:/ /www- 

I sk. icrr . u-tokyo . ac . j p/ doc/ sk/pu b 



[28] K. Inoue, in TAUP '97, Proceedings of the 5th International Workshop on Topics in 
AstroParticle and Underground Physics (Gran Sasso, Italy, 1997), to appear in Nucl. 
Phys. B (Proc. Suppl.); also available at the same URL as in ref. ||27|| . 
Y. Suzuki, in Neutrinos in Astro, Particle and Nuclear Physics, Proceedings of the 19th 
International School of Nuclear Physics (Erice, Italy, 1997), to appear; also available at 



the same URL as ref. 27 



[30] Y. Itow, in the Proceedings of the Topical Conference "Physics of Leptons," SLAG 
Summer Institute 1997; (Stanford, GA, 1997), to appear; also available at the same 
URL as ref. 0. 

[31] H. Sobel, in the Proceedings of the Aspen Winter Conference on Particle Physics: Fron- 
tiers in Particle Physics (Aspen, GO, 1998), to appear. 
[32] I. Pomeranchuk, as cited in Gribov and Pontecorvo [0. 

[33] G. L. Fogh, E. Lisi, D. Montanino, and G. Scioscia, Phys. Rev. D 56, 4365 (1997). 
[34] J. N. Bahcall, M. Kamionkowsky, and A. Sirhn, Phys. Rev. D 51, 6146 (1995). 



14 



FIGURES 



FIG. 1. Vacuum oscillation solutions to the solar neutrino deficit in the usual mass-mixing 
plane, as derived from a fit to the data of Table I. The four regions A, B, C, and D, are allowed at 
95% C.L. 

FIG. 2. Energy spcctrTim of recoil electrons in Borexino. The main components are shown 
separately. The smearing effect of the energy resolution is also shown. The arrows indicate the 
prospective energy window assumed in the analysis. 

FIG. 3. SuperKamiokande (SK) experiment. Map of the four vacuum oscillation solutions A, 

B, C, and D in the plane charted by A(T) / (T) (the fractional deviation of the mean electron kinetic 
energy) and by /i — e (the deviation of the first Fourier component from its standard value) . Notice 
the strong correlation between these two variables, which is induced by the L/E dependence of the 
oscillation probability. The star at the origin corresponds to the standard (no oscillation) case. The 
SuperKamiokande datum on A{T)/{T) is also shown; it disfavors solution C at > 4(7 and solution 
D at > 6a. The part of the solution A which is favored by the SK datum at la corresponds to an 
expected Fourier amplitude /i — e ~ 0.5-1%. 

FIG. 4. Borexino (BX) experiment. Map of the four vacuum oscillation solutions A, B, C, and 
D, in the plane charted by the first two Fourier components. The star at the origin corresponds to 
the no oscillation case. There is some overlap among the four solutions. 

FIG. 5. SuperKamiokande (SK) vs Borexino (BX). Map of the four vacuum oscillation solutions 
A, B, C, and D in the planes charted by the SK observables A{T)/{T) and f\ — e (ordinates) and 

by the BX observables fi — £ and /2 (abscissae). Future measurements of the /„'s in SK and 
BX, together with increasingly accurate A{T)/{T) data from SK, are expected to spot one of the 
solutions. 
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VACUUM SOLUTIONS 
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FIG. 1. Vacuum oscillation solutions to the solar neutrino deficit in the usual mass-mixing 
plane, as derived from a fit to the data of Table I. The four regions A, B, C, and D, are allowed at 
95% C.L. 
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Borexino, standard spectrum 
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FIG. 2. Energy spectrum of recoil electrons in Borexino. The main components are shown 
separately. The smearing effect of the energy resolution is also shown. The arrows indicate the 
prospective energy window assumed in the analysis. 
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FIG. 3. SuperKamiokande (SK) experiment. Map of the four vacuum oscillation solutions A, B, 
C, and D in the plane charted by A{T) / (T) (the fractional deviation of the mean electron kinetic 
energy) and by /i — e (the deviation of the first Fourier component from its standard value) . Notice 
the strong correlation between these two variables, which is induced by the L/E dependence of the 
oscillation probability. The star at the origin corresponds to the standard (no oscillation) case. The 
SuperKamiokande datum on A{T)/{T) is also shown; it disfavors solution C at > 4(T and solution 
D at > 6cr. The part of the solution A which is favored by the SK datum at la corresponds to an 
expected Fourier amplitude /i — e ~ 0.5-1%. 
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FIG. 4. Borexino (BX) experiment. Map of the four vacuum oscillation solutions A, B, C, and 
D, in the plane charted by the first two Fourier components. The star at the origin corresponds to 
the no oscillation case. There is some overlap among the four solutions. 
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SUPERKAMIOKANDE vs BOREXINO 




FIG. 5. SuperKamiokande (SK) vs Borexino (BX). Map of the four vacuum oscillation solutions 
A, B, C, and D in the planes charted by the SK observables A{T) / (T) and /i — e (ordinates) and 
by the BX observables fi — e and /2 (abscissae). Future measurements of the /n's in SK and 
BX, together with increasingly accurate A{T)/{T) data from SK, are expected to spot one of the 
solutions. 
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